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Abstract 

We compute long distance effects on the photon energy spectrum in inclusive 
radiative decays of B mesons using light-cone expansion and heavy quark ef- 
fective theory. We show that for sufficiently high photon energy the leading 
nonperturbative QCD contribution is attributed to the distribution function. 
The distribution function is found to be universal in the sense that the same 
distribution function also encodes the leading nonperturbative contribution to 
inclusive semileptonic decays of B mesons at large momentum transfer. Some 
basic properties of the distribution function are deduced in QCD. Ways of 
extracting the distribution function directly from experiment and their impli- 
cations are discussed. The theoretically clean methods for the determination 
of \Vts\ are described. 
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I. INTRODUCTION 



The study of the inclusive radiative decay B Xg'y, where Xg is any possible final 
state of total strangeness —1, is a broad subject with many areas of investigation. It is a 
fiavour changing neutral current process, which is forbidden at tree level, only proceeding 
through what is called a electroweak penguin diagram in the Standard Model. This is 
a one-loop graph inducing b s+photon with a W boson and a quark, predominantly the 
top quark, in the loop. Measurements of this rare process provide one of the most stringent 
experimental tests of the Standard Model at one-loop level. It can be used to determine the 
fundamental Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix element \Vts\ and test its 
unitarity. It may offer new insight into the nature of confinement. 

Meanwhile, because the decay is highly suppressed in the Standard Model, it could 
be particularly sensitive to new physics beyond the Standard Model. Measurements of 
B — > Xgj decays would impose constraints on new physics models. In particular, new con- 
tributions from nonstandard particles replacing the standard model particles in the loop 
can be detected prior to a direct production of such new particles at much higher ener- 
gies. A clear deviation from standard model expectations would signal new physics from 
super symmetry, charged Higgs scalars, anomalous WW'j couplings, etc. 

Observation of the inclusive radiative decay B X^'j was reported in 1995 by the CLEO 
Collaboration, obtaining the first measurement of the branching ratio of (2.32±0.57±0.35) x 
10~^ 0- The ALEPH Collaboration has reported a measurement of the corresponding 
branching ratio of beauty hadrons produced at the Z resonance, obtaining the branching 
ratio of (3.11 ± 0.80 ± 0.72) x 10"^ 13 . Recently, the CLEO Collaboration has improved 
their measurement using 60% additional data and improved analysis techniques, yielding 
a new preliminary result of the branching ratio of (3.15 ± 0.35 ± 0.32 ± 0.26) x 10~^ [lll| . 



Copious data samples from the rare decays B Xg'y will soon be collected at B factories. 
The accuracy of the B Xg'y measurement will be significantly improved shortly. 

Significant progress has been made in the theoretical prediction of the branching ratio 
for B — > Xs'y in the Standard Model. The next-to-leading order (NLO) calculation of the 
perturbative QCD corrections has been completed recently, combining the calculations of 



the matching conditions of the Wilson coefficients |T2-Q5i, the matrix elements [16-18], and 



the anomalous dimensions iTPIl. This achievement leads to a substantial reduction of the 



large renormalization scale dependence of the leading order result |20-28|. The leading 



power corrections in 1/ml |^,30| and 1/m^ |31-05[] have been calculated in the context of 



the heavy quark expansion. In addition, the leading electroweak radiative corrections have 
recently been investigated []36| , ^ |. 

It is, however, not sufficient to have a reliable calculation of the total decay rate. Ex- 
perimentally, in order to suppress backgrounds from other S-decay processes, only events 
in the high energy region of the photon energy spectrum in 5 — Xg'y decays have been 
used to measure the branching ratio. An accurate and reliable theoretical description of the 
photon energy spectrum is essential in order to perform a fit to the experimental data and 
extrapolate to the total decay rate. The study of the photon energy spectrum is also quite 
interesting in its own right, since it would probe the decay dynamics in more detail than 
the total rate. Strong interaction effects can significantly modify the photon energy spec- 
trum. The following fact well illustrates the situation. It is well known that the quark-level 
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and virtual gluon exchange processes 6 — 57 generate a trivial photon energy spectrum 
— a discrete line at = mfe/2 in the 6-quark rest frame. It is two distinct effects, gluon 
bremsstrahlung and hadronic bound state effects, that spread out the spectrum. Gluon 
bremsstrahlung results in a long tail in the photon spectrum below mi,/ 2. Bound state ef- 
fects leads to the extension of phase space from the parton level to the hadron level, also 
stretches the spectrum downward below mf,/2 and is solely responsible for populating the 
spectrum upward in the gap between the parton level endpoint = mb/2 and the hadron 
level endpoint E^ = Mb/2. 

The perturbative QCD corrections to the photon energy spectrum are known at order 
as [0 and the order a^Po contribution is computed very recently in the perturba- 
tion expansion of the matrix elements. The relevant Wilson coefficients in the effective 
weak Hamiltonian are known to next-to-leading logarithmic accuracy |T^ in renormaliza- 
tion group improved perturbative QCD. Currently nonperturbative QCD effects comprise 
the potentially most serious source of theoretical error in the photon energy spectrum. Long 
distance QCD effects are relevant for improving the accuracy of an analysis and understand- 
ing theoretical uncertainties in the Standard Model, as well as for searching for new physics 
by the confrontation of the standard model predictions with the data. In view of this, it is 
important to address the issue as rigorously and completely as possible. The Fermi motion 
of the b quark in the B meson has been taken into account in WB^ by using the phenomeno- 



logical model by Altarelli et al. A more fundamental treatment of bound state effects 



based on the heavy quark expansion has been developed by resumming an infinite set of 
leading- twist operators into a shape function |]40|-^ . Recently, an analysis of long distance 
effects on the photon spectrum has been performed along these lines, including for the 
first time the full NLO perturbative QCD corrections. 

In this paper we will compute the long distance QCD contributions to the photon energy 
spectrum using light-cone expansion and heavy quark effective theory (HQET). These ideas 
and techniques have originally been exploited in the theoretical description of inclusive 
semileptonic decays of beauty hadrons |^3[. This approach has been related |^ to the 



above-mentioned similar approach based on the resummation of the heavy quark expansion 
In this paper we extend our previous analyses to inclusive radiative decays of B 



mesons. We strive to present our approach in a more detailed and complete form. 

The fact that beauty hadrons are the heaviest hadrons actually formed since the top 
quark is too heavy to build hadrons confers a special role to beauty hadrons. The beauty 
hadron decays involve two large scales: the heavy beauty hadron mass at the hadron level 
and the heavy b quark mass at the parton level, which are much greater than the energy 
scale Aqcd which characterizes the strong interactions. The two large scales give rise to 
two techniques for dealing with nonperturbative QCD. Because of the heaviness of the 
decaying hadron, the decay dynamics may be dominated by the space-time separations in 
the neighborhood of the light cone. The light-cone expansion provides a formal and powerful 
way of organizing the nonperturbative QCD effects and singling out the leading term. On 
the other hand, the heavy quark mass provides a large limit to construct an effective theory 
describing the heavy quark interacting with the gluons in the heavy hadron. This so-called 
heavy quark effective theory has new (approximate) symmetries that were not manifested 
in the full QCD Lagrangian and sets another framework for organizing and parametrizing 
nonperturbative effects, which relates various phenomena (e.g., the hadron spectroscopy and 
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weak decays of hadrons containing a single heavy quark) to a common set of parameters, 
so that it has great predictive power. It is certainly useful to calculate long distance QCD 
effects on the photon energy spectrum in B ^ X^'y decays using light-cone expansion and 
heavy quark effective theory. 

The paper is organized as follows. In Sec. 11, we generalize the results of Ref. |^3[ and 
construct the light-cone expansion for the photon energy spectrum in B —>■ Xgj decays. We 
show that for sufficiently high photon energy, the leading nonperturbative QCD contribu- 
tion is attributed to a distribution function. We define a new, gauge-invariant distribution 
function. The connection between the nonperturbative QCD effects in radiative and semilep- 
tonic inclusive B decays is subsequently established in Sec. 111. For the sake of completeness, 
we study the properties of the new distribution function in Sec. IV. We discuss the direct 
extraction of the distribution function from experiment and precise determinations of the 
CKM matrix elements in Sec. V. We hope to motivate the measurement of the distribution 
function by pointing out its importance in inclusive B decays and the ways to measure it. 
Section VI briefly summarizes the main results. 



II. LIGHT-CONE EXPANSION FOR B X^j DECAYS 

We shall study inclusive radiative decays B{P) Xs{px)'y{P'y)- A suitable framework 
to do that is an effective low-energy theory, obtained by integrating out the heavy particles 



which in the Standard Model are the top quark and the W boson Since we shall concen- 
trate on long distance effects, it suffices to work at lowest order in electroweak interactions 
ignoring for the moment perturbative QCD and electroweak radiative corrections. To lowest 
order B X^'-f decays are governed by the effective weak Hamiltonian involving a magnetic 
dipole operator: 



where the coupling constant 



Kfr = nsaf^^RbF^,, (2.1) 



gauges the strength of the b ^ S'j transition, R = {1 + 75)/2 is the right-handed projection 
operator, and F^^^, is the electromagnetic field strength tensor. In Eq. (p.2|), Vtb and Vts are 
the CKM matrix elements, the function CY'\My\/) depends on the masses of the internal 
top quark and W boson, and takes the form P 

^^^^ - 4(x,-l)4^^'^* ~ 24(a:, - 1)3 ^^'^^ 

with Xf = rri^/M^. We use to denote the b quark mass and Mb to denote the B meson 
mass. It can be shown that contributions of the strange quark mass to the decay rate are 
suppressed by the very small factor m^/mf. Hence we neglect the strange quark mass unless 
noted otherwise. 

The decay rate for B Xg'j is given by 
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d^ = 7^ S^,:j, E(2^)'^'(-P-P7-Px)|(X.(px)7(p.,6)|Kff|i?(P))P, (2.4) 

where we sum over all possible final states of total strangeness —1 as well as the two 
transverse polarizations of the photon. We adopt the standard covariant normalization 
{B\B) = 2£'5(27r)^(5^(0) for the B meson state. The decay rate can be expressed in terms 
of a current commutator taken between the B meson states, which incorporates all nonper- 
turbative QCD physics of the weak radiative i?-meson decay. We find 

with the generalized current J^(y) = b{y)['y^, ^^]Ls{y) , where L = (1 — 75)/2 is the left- 
handed projection operator. 

The commutator in Eq. (|2.5| ) has to vanish for space-like separations < due to 
causality. Moreover, the behavior of the integral in Eq. (|2.5| ) is determined by the integrand 
in domains with less rapid oscillations, i.e., ~ 1, which implies y'^ ^ l/-^'^- Combining 

these results we find that the dominant contribution to the integral in Eq. ( p.5|) results from 
the space-time region < y^ ^ l/-^7- This implies that for sufficiently high photon energy, 
^ Aqcd; the space-time separations in the neighborhood of the light-cone = 
dominate the decay dynamics. It is then appropriate to construct a light-cone expansion to 
calculate the photon energy spectrum in the high energy region. 

The decay produced strange quark propagating in a background gluon field obeys the 
anticommutation relation 

{s{x),s{y)} = ms)i^s{x -y)U{x,y), (2.6) 

with the Wilson link 

U{x, y) = VeMWs r dz^'A^iz)] (2.7) 
Jy 

between the quark fields at y and x, where V denotes path ordering. Here the strange quark 
mass nis 7^ is retained. The background gluon field is that generated by the remnants 
of the decaying B meson. Aq{y) is the Pauli- Jordan function of the form [^] 



A,iy) = / d^ke-^'^ yeik^ik' - mj) 



= -7r^{y')S{y') + -^e{yy{y')Um,^^). (2.8) 
27r 47rvy 

Here Ji{z) is the Bessel function of order 1. From Eq. ( |2.8| ) we see that Ag(?/) is singular 
on the light cone and would also select out light-cone contributions were it not for high 
frequency variations in the phase. Using the anticommutation relation (|2.6| ), we find 



[Jliy). J'm = iQp,^p,,[id-iAs{y)]mi^RU{^.y)h{y). (2.9) 

Substituting Eq. (OI) in Eq. (O), it follows that 
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= ^ (2^)3^2^ ^^""^"^ / ^'2/^*'"'[^^°^^^(^)K^IK0)7^t/(0,y)6(y)|5). (2.10) 

Note that the matrix element of b{0)-y'^j5U{0,y)b{y) between the B meson states vanishes 
by parity invariance in the strong interactions. 

Now let us consider the matrix element of the bilocal operator in Eq. ( p.lOl) . The general 
tensor decomposition of it reads 

{B\b{Q)^^UiO,y)b{y)\B) = 2[P^F{y\y P) + y^G{y\y ■ P)], (2.11) 

where F{y'^, y ■ P) and G{y'^, y ■ P) are in general functions of the two independent Lorentz 
scalars, y'^ and y ■ P. Contracting both sides of Eq. ( f^.llj ) with yp gives 

(5|6(0)^t/(0, y)b{y)\B) = 2[y PF{y\ y ■ P) + y^G{y\ y ■ P)]. (2.12) 

In the domain of interest where \P'y ■ y\ ~ 1 and < ^ 1/i?^, we infer from Eq. (p.l2|) 
that the contribution of G{y'^,y ■ P) is suppressed by a factor for high photon energy 
considered here, compared to that of F{y'^,y ■ P). In addition, we can make a light-cone 
expansion for F{y'^,y ■ P), which is justified for high photon energy as discussed above. 



n=0 



d-F{y',y.P) _ ^2.13) 

y2=o 



dy 



2n 



The rath term in the light-cone expansion is suppressed by XjE'^ . Therefore, in the leading- 
twist approximation we then have 

(5|K0)7''f/(0,y)%)|i?) = IpfFiy" = 0,2/ ■ P). (2.14) 

The contributions from F{%j^ 0,y ■ P) and G{y'^, y ■ P) represent higher twist effects. They 
become important only at the subleading level, a question that will not be dwelt upon here. 
On the light cone P(?/^ = 0,y ■ P) can be projected out of the general decomposition 

(my- 

Fiy' = O,y.P) = ^^{Bm0iO,y)biy)\B)\y.=,. (2.15) 
The Fourier transform of P(?/^ = 0,y ■ P) defines the distribution function 



fiO = l^ fdiyP)e^^y''F{y' = 0,y.P) 
Zn J 



= i/ ^^r^e^^'-''(^IKO)^f/(0,y)6(y)|P)|,.=o- (2.16) 

The Wilson link is associated with the background gluon field, which ensures gauge invari- 
ance of the distribution function. The leading nonperturbative QCD contribution contained 
in P(?/^ = 0,y ■ P) is equivalently translated into the distribution function. 
Using the inverse Fourier transform 
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F{y^ = 0,y.p)=J d^e-'^y^'fiO , (2.17) 

substituting Eqs. ( p.8|) and ( p.l4|) in Eq. (|2.10|) , and then carrying out the phase space 
integration, we arrive at the photon energy spectrum in the B rest frame 

^^(^ ^ ^-^^ = ^^iv,,v,:ncf (M^)pE^/f^V (2.18) 



Equation ( p^.l8| ) holds in the leading twist approximation, which is expected to be a good 



approximation for the energy of the emitted photon of around 1 GeV and above. The precise 
shape of the spectrum near the lower endpoint i?^ = is not available to us, where the light- 
cone expansion is not applicable and higher twist terms give contributions of the same order 
as the leading twist term. However, given the experimental P- |ll|] and theoretical ||T6| , |42|j37 



indications that the spectrum for < 1 GeV appears to be vanishingly small, a more 
accurate account of the rather small overall nonperturbative contributions in this region 
seems numerically unimportant. This is not unexpected since the photon energy spectrum 
stemming from the quark-level and virtual gluon exchange processes would only concentrate 
at E^ = mh/2 ^ 2.45 GeV and gluon bremsstrahlung and hadronic bound state effects smear 
the spectrum about this point, but most of the decay rate remains at large values of E^. 
For practical purposes, this approach would nonetheless provide a realistic description of the 
photon spectrum over the full kinematic range < E^ < Mb/2 = 2.64 GeV. Of course, one 
should keep in mind that for the photon spectrum near the lower endpoint one has nothing 
from the approach to be verified or falsified. 



Finally, let us discuss the relation of our approach to the approach advocated in [4C-42 
based on the resummation of the heavy quark expansion. By assuming E^ = mi,/ 2 for the 
factor E^ in Eq. (|2.18|) we can reproduce their formulas for the B photon energy 



spectrum obtained in pT| , ^ . In that case, the photon energy is fixed to the value in the 
free quark limit, instead of varying in its entire kinematic range from to Mb/ 2. 



III. UNIVERSALITY 

In this section we compare the inclusive radiative decay with the inclusive semileptonic 
decays of B mesons B Xflv (/ = m, c) in order to decipher the universal, process- 
independent structure in the relevant realm of hadron physics. 

Inclusive semileptonic decays of B mesons at large momentum transfer carried by the W 
boson are also governed by the light-cone dynamics 
on the decays reside in the hadronic tensor 



The nonperturbative QCD effects 



W. 



1 

2^ 



d'ye^'^-y{B\[Uy),jliO)]\B), 



(3.1) 



where the charged weak current j^{y) = /(?/)7^(l — 75)6(1/). Calculating the charged weak 
current commutator in the same way as in the last section yields 



{B\ \j,{y),jl{0)] \B) = 2{S,^,^ - ie,^,p) [9" A^(y)] (i?|6(0)7^f/(0, y)h{y)\B) 



(3.2) 
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where S^aup = gtiagiyis + g^l(}gua — giiiygaf3- We observe here the very same matrix element, 
{B\b{0)'y^U{0,y)b{y)\B), entering the decay rate for the inclusive semileptonic decays of B 
mesons. 

The twist expansion of the bilocal operator matrix element then leads to the expression 
of the differential decay rate in the B rest frame in terms of the distribution function /(^) 
defined in Eq. ( |2.16|) : 

"^^'f-f^-' ^ {fME^M. - f) - iU - ^-)} . (3.3) 

where 



go ± A/lqP + ml 



B 



and mj is the mass of the decay produced quark, which is the up (charm) quark with 
/ = u{c). All lepton masses have been neglected. Therefore, the distribution function 
describes leading long-distance QCD effects not only in inclusive radiative B decays, but 
also in inclusive semileptonic B decays, and is thus a universal function. This universality 
originates from the fact that the primary object of analysis in long distance effects is the 
same bilocal operator matrix element dictated by the light-cone dynamics. 

The b quark distribution function /(^) defined in Eq. ( p.l6| ) differs from the definition 
given in . The difference arises as the distribution function defined here is gauge invariant 



and its inverse Fourier transform is exactly equal to the leading twist term -F(y^ = 0,y ■ P), 
whereas the distribution function of is defined in the light-cone gauge and its inverse 



Fourier transform is equal to F[y'^ = 0, y-P) in the leading twist approximation. However, we 
note that the expressions for the inclusive semileptonic B decay rates derived previously |^ 



in terms of the distribution function remain formally the same, if the distribution function 
defined in this paper is used. This is exemplified by Eq. (|3.3|) for the triple differential decay 
rate. 



IV. PROPERTIES OF THE DISTRIBUTION FUNCTION 

The b quark distribution function is an important physical quantity which summarizes 
leading long distance effects on inclusive B decay processes. Equation ( ^.18| ) shows that the 
photon energy spectrum depends strongly on the distribution function. In Sec. II we defined 
the b quark distribution function in QCD by modifying the definition given in Ref. In 



this section we re-derive several important properties of the distribution function. Since it 
is gauge invariant, it can be calculated in any gauge. For simplicity, we shall choose the 
light-cone gauge in the following discussion so that the Wilson link becomes the identity 
operator. 

It is helpful to introduce the null vector = (1,0,0,-1) at this step, such that the 
light-like vector y^^ = tn'^ with t being a parameter. Using this notation, the distribution 
function defined in Eq. ( p. 161 ) can be rewritten as follows 
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/(O = ^J dte^^'--''{B\bU0)b4tn)\B), (4.1) 

where 6+ = is the "good" component projected out of the b quark field, with the 
projection operator P+ = (1 + 7*^7^)72. Inserting a complete set of hadronic states between 
quark fields and translating the tn dependence out of the field, we find 

fiO =Y.^l^-iP-^P- Pm)]Km|6+(0)|S)p. (4.2) 

m 

So the distribution function obeys positivity. The state |m) is physical and must have 
< < Eb, thus /(^) = for ^ < or ^ > 1. Therefore, the support of the distribution 
function reads < ^ < 1. Moreover, one can observe from Eq. (f4.2|) that /(^) is the 
probability of finding a b quark with momentum inside the B meson. This is the familiar 
probabilistic interpretation of the parton model for inclusive B decays. 



From Eqs. ( p.l7|) and (|2.11| ) it is straightforward to show that 

f'dUiO = Fiy = 0) = l, (4.3) 
Jo 

because b quantum number conservation implies 

{B\lrf^'b\B) = 2P^. (4.4) 

Thus the distribution function is exactly normalized to unity, which does not get renormal- 
ized as a consequence of b quantum number conservation. 

To understand B meson bound state effects, it is instructive to consider first of all the 
form of the b quark distribution function in the free quark limit. In this limit, the B meson 
and the b quark in it move together with the same velocity: pt/tTif, = P/Mb = v and the 
free Dirac field b{y) = e~*^'^''6(0), so from Eq. ( p^.l6D it follows that 

hee{0=Si^-mb/MB). (4.5) 

Substituting Eq. ( [4.5| ) in Eq. ( p.l8| ), it consistently reduces to the free quark decay spectrum 
in the rest frame of the b quark 

dT{b -> s-f) Glani ^ nib , . 

dE, = ^2;^'^*^^-' '^^ ^^^^ - T^' (^-^^ 

which is a discrete line at E^ = mi,/ 2. Reversely, the physical photon spectrum is ob- 
tained from a convolution of the hard perturbative spectrum with the soft nonperturbative 
distribution function /(O- 

aL^ Jo dE^ 

The bound state smearing of the photon energy spectrum is then refiected in the deviation 
of the true distribution function from the delta function form. The quantitative analysis of 
this deviation is the subject of the rest of this section. 
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Additional properties of the distribution function can be deduced by means of operator 
product expansion and heavy quark effective theory. The derivation proceeds in essentially 
the same manner as in Ref. but for the different definition of the distribution func- 
tion. We present the corresponding derivation for the new distribution function below for 
completeness. 

Since the b quark inside the B meson behaves as almost free due to its large mass, relative 
to which its binding to the light constituents is weak, one can extract the large space-time 
dependence 

6(y)=e-^-^-^6,(l/). (4.8) 

A Taylor expansion of the field in a gauge-covariant form relates the bilocal and local oper- 
ators. This leads to an operator product expansion 



mi^Ky) = e-'"^'""" E ^^A^i ■ ■ ■ yM^)i^k^'' ■ ■ ■ k''-\{o), (4.9) 



n=0 



where fc^ = iDfj_ = i{d^ — igsA^) and the symbol {■ ■ ■} means symmetrization with respect 
to the enclosed indices. Using Lorentz covariance one can express the matrix element of the 
local operator on the right-hand side of Eq. ( |4.9| ) between the B meson states in terms of 
the B meson momentum: 

(B|6^(0)7^A;^^i ■ ■ ■ P">fe^(0)|5) = 

n 

2{CnQP'^P^' . . . + ^ M^Cnig'^'^'P^' ■ ■ ■ P^-ip^»+i . . . P^") 

i=l 

+ terms with g^^'^^3 . (4.10) 

The terms with g^^^^ drop out on the light cone. Substituting Eqs. ( [4.9D and ( [4.10[ ) in 
Eq. ( ^161) yields 

/(0 = i:^^c7>")(e-^). (4.11) 

Therefore we obtain the moment relation for the distribution function 



Mn{mb/MB) = CnO, (4.12) 

where the nth-moment about a point ^ of the distribution function is in general defined by 

MniO= f'da^-OViO- (4.13) 
Jo 

By definition, Mo(^) = Cqq = 1. The moment relation ( |4.12| ) relates the moments of the 
distribution function to the matrix elements of the local operators in Eq. (|4.10| ). 

We invoke heavy quark effective theory to evaluate further expansion coefficients in 
Eq.( [4.10|) . In this effective theory the QCD 6-quark field b{y) is related to its HQET coun- 
terpart h{y) by means of an expansion in powers of l/mf,: 
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1 + 



2mb 



ml 



Ky). 



The effective Lagrangian takes the form 



-HQET 



hiv ■ Dh 



2mh 



mt 



(4.14) 



(4.15) 



where QsG"^^ = ^[-D", -D^] is the gluon field-strength tensor. Only the first term in Eq. (f4.15|) 
remains in the rrify oo limit, which respects the heavy quark spin-flavor symmetry. The 
other two terms give the l/rrib corrections: the second term violates the heavy flavor sym- 
metry, while the third term violates both the spin and flavor symmetries. Using the method 
described in Refs. [^7|-^ to relate matrix elements of local operators in full QCD to those 
in the HQET, the expansion coefficients Cni in Eq. (|4.10|) can be expressed in terms of the 
HQET parameters. The nonperturbative QCD effects can, in principle, be calculated in a 
systematic manner. In this formalism the moment Mn{fnb/ Mb) is expected to be of order 
(AqcD/w^fc)"- A few coefficients have been evaluated and the results are |^ 



Cio 

C20 
C21 



3Mb 



3Ml 



ir, + 0(A^CD/<), 



G 



22 



0, 



(4.16) 
(4.17) 

(4.18) 
(4.19) 



where Eh = Kh + Gb and Kb and Gb are the dimensionless HQET parameters of order 
(AqcD/w^b)^, which are often referred to by the alternate names Ai = —2mlKb and A2 = 
—2mlGb/3, defined as 



Ai 
A2 



1 



2Mb 
1 



12M 



{B\h{iDfh\B), 
-{B\hg,Go.(3a''^h\B). 



B 



(4.20) 
(4.21) 

(M|,-M|)/4 = 



The parameter A2 can be extracted from the B* — B mass splitting: A2 
0.12 GeV^. The parameter Ai suffers from large uncertainty. 

Thus two sum rules for the distribution function can be derived according to the moment 
relation ([4.12|) . They determine up to order (AQCD/''^b)^ the mean value fi and the variance 
(T^ of the distribution function, which characterize the location of the "center of mass" of 
the distribution function and the square of its width, respectively: 



rUb 



M, 



B 



a 



Mr 



(- 


5Eb\ 




3 j' 


2Kb 


-n 


3 





(4.22) 



(4.23) 
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with the definitions 



/i = Mi(0)=|+Mi(O, (4.24) 
= M2(/x) = M^ii) - M^iO. (4.25) 

The mean value and variance specify the primary shape of the distribution function. There- 
fore, our evaluation comes to the conclusion that the distribution function /(^) is sharply 
peaked around ^ = ~ mb/Ms close to 1 and its width of order Aqqd/Mb is narrow. 

The sum rules given in Eqs. (|]22D and (^1231) quantify the deviat ion of the true dis- 



tribution function from the delta function form due to bound state effects. Choosing the 
parameters rrib = 4.9 GeV and Ai = —0.5 GeV^ for the purpose of orientation, we ob- 
tain from Eqs. ( |4.22| ) and ( ^.23| ) fi = 0.933 and cx^ = 0.006. By contrast, the mean value 



H = rrih/MB = 0.928 and the variance = in the free quark limit. 

The results for the mean value and the variance given in Eqs. ( |4.22| ) and ( [4.23] ) are at 



variance with those of due to the different definitions of the distribution function. How- 
ever, the corresponding numerical shifts are found to be so small that the phenomenological 
impacts of these differences are negligible. 

Nonperturbative QCD methods such as lattice simulation and QCD sum rules could help 
determine further the form of the distribution function. The distribution function could also 
be extracted directly from experiment, as we shall discuss below. 



V. MEASUREMENTS OF THE DISTRIBUTION FUNCTION AND 
DETERMINATIONS OF THE CKM MATRIX ELEMENTS 

The universality of the distribution function discussed in Sec. Ill enhances the predictive 
power of the approach: the distribution function can be extracted from measurements of 
one process and then used to make predictions in all other processes in a model-independent 
manner. In this section we start by exploring how to extract the distribution function from 
experiment. Then we will investigate the theoretically clean methods for the determination 
of the CKM matrix element \Vts\- 

The spectrum result, Eq. (|2.18| ), can be cast in the following form: 

,,.,2f(^EA _ 1 dT{B -> X,7) 



\MbI Glaml\cf\Mw)?E^^ dE, 



This immediately implies that the distribution function can be extracted directly from a 
measurement of the photon energy spectrum upon implementing perturbative QCD and 
electroweak radiative corrections, which have so far been ignored in this paper. The exper- 
imental data on the photon energy spectrum are already available [p|-pn|, but limited by 



statistics for a meaningful extraction of the distribution function. Forthcoming very large 
data samples from high-luminosity B factories promise a direct extraction of the distribution 
function with reasonable accuracy. Such a determination of the distribution function would 
help substantially to reduce the theoretical uncertainties in the descriptions of both semilep- 
tonic and radiative inclusive decays of B mesons. The distribution function extracted from 
B — > can be applied, for example, in the calculations of the lepton energy spectrum and 
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the hadronic invariant mass spectrum in charmless inclusive semileptonic decays B Xuijy, 
so that the precision of the \Vub\ determination from these spectra can be improved. 

The distribution function can also be extracted directly from the measurements of the 
differential decay rates as a function of the scaling variable ^+ [note that ^+ is different 
kinematic variable for B —>■ Xulv and B —>■ Xciiy, defined in Eq. (|3.4| )] in the inclusive 
semileptonic decays of B mesons [Q: 



where r/ = mj/MB (/ = u,c) and <l>(x) = 1 — 8x^ + 8x^ — — 24x^lnx. Such a deter- 
mination of the distribution function would also benefit the theoretical descriptions of both 
semileptonic and radiative inclusive decays of B mesons. In particular, once such an ex- 
traction of the distribution function is available, it can be applied as an independent input 
in turn in the analysis of the B Xg'y photon energy spectrum. This would substantially 
improve the standard model predictions and increase the sensitivity to new physics. There- 
fore, this extraction of the distribution function is valuable in the search for additional decay 
mechanisms beyond the Standard Model in B ^ Xg'j decays. The experimental technique 
of neutrino reconstruction could well make this way of extracting the distribution function 
experimentally feasible. If the neutrino can be reconstructed kinematically by inferring its 
four- momentum from the missing energy and missing momentum in each event, then it is 
possible to measure the scaling variable 

It is important to note that the B Xs'j photon energy spectrum and the B Xfiu 
spectra dr/dC,+ share a common feature. Namely, they offer the intrinsically most sensitive 
probe of long-distance strong interactions because these spectra correspond to a discrete 
line solely on kinematic grounds in the absence of gluon bremsstrahlung and long-distance 
strong interactions. Indeed, our calculation based on the light-cone expansion shows that 
they are explicitly proportional to the nonperturbative distribution function. Therefore, the 
shapes of these spectra directly reflect the inner long-distance dynamics of the reactions and 
measurements of these spectra are idealy suitable for direct extraction of the distribution 
function from experiment. This salient feature also makes the 6 — > u scaling variable ^+ 
unique to give a very efficient discrimination between B X„£z/ and B — > Xciiy events ||50| , 
even better than the hadronic invariant mass. 

Given the difficulty in distinguishing t ^ s decays from the dominant t —>■ b decay mode, 
it is hard to determine \Vts\ from studies of top quark decays. Based on the results obtained 
in this paper and in Ref. |50|, we propose a new strategy to extract \Vts\ which is largely free 
of hadronic uncertainties. The idea is to use the known normalization of the distribution 
function or the cancellation of the distribution function in the ratio of the decay rates to 
eliminate the dominant hadronic uncertainties. The most straightforward and best way to 
eliminate the dependence on the distribution function is to resort to, again, the B — > X^^ 
photon energy spectrum and the B — > X^^c^v spectra dT/d^j^ since they are proportional to 
the distribution function. 

Integrating Eq. (|5.1|) over and using the normalization condition ( [4.3|) yields 

r*i2 _ 47r4 /"^^^/^jn. I dT{B ^ Xs-f) 



\VtbV:y = TTT. / dE^— ' ,^ — (5.3) 
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Thus the known normahzation of the distribution function allows an almost model indepen- 
dent determination of |VJbVJs| from a measurement of the weighted integral of the photon 
energy spectrum. This determination can be taken as a measurement of \Vts\ by using 
\Vtb\ = 0.9993 from unitarity of the CKM matrix, which holds to a very high accuracy []5l 



The advantage of this determination of \Vts\ is that the dominant hadronic uncertainty has 
been avoided, which may provide one of the most precise determinations of \Vts\- 

As in the case of 5 Xg'j discussed above, one can get rid of the distribution function 
by integrating Eq. ( p.2|) over ^+ |[50[] : 



fb\ 



1927r3 

Gp4 



dT{B ^ Xfiu) 



(5.4) 



Thus, the precise determinations of \Vub\ and |V"cfe| can be obtained from the measurements 
of the weighted integrals of the differential decay rates as functions of ^+ in i? ^ X^ii^ and 
B Xciiy, respectively. 

Alternatively, when we take the ratio of the differential decay rates, the distribution 

and 



function cancels. From Eqs. (|5Tl 



we obtain 



fb 



ttMb 



3aml\cP{Mw)\ 



dT{B Xs-f)/dE^ 



dT{B ^ Xf£u)/d^^ 



(5.5) 



E^=Mb^+/2 



which may be useful to provide a theoretically clean determination of |Vjs/V/b|. By the same 
token, an analogous expression has been derived in for B Xu^c^u, which can be used 
to measure |V^fe/K;fe| to a high precision. These determinations of the CKM matrix elements 
rely on the universality of the distribution function, in contrast to the method described 
above by virtue of the known normalization of the distribution function. The compatibility 
between the various measurements will test the universality of the distribution function in 
the inclusive semileptonic and radiative decays of B mesons within the Standard Model. 



VI. SUMMARY 

In this paper we have calculated the long distance effects on the photon energy spectrum 
in B ^ Xg'y decays. We have demonstrated on the basis of light-cone expansion that the 
leading nonperturbative QCD contribution in inclusive radiative B decays with emission 
of a sufficiently high energy photon resides in the distribution function. The distribution 
function is defined by Fourier transformation of the matrix element of the non-local b quark 
operators separated along the light cone. We have found that the distribution function 
is universal in the sense that the same distribution function also summarizes the leading 
nonperturbative QCD contribution in inclusive semileptonic B decays at large momentum 
transfer. 

Although long-distance strong interactions responsible for the distribution function pre- 
clude a complete calculation of it at present, we have deduced some of its basic properties 
in QCD. The distribution function is gauge invariant and obeys positivity. It has a support 
between and 1 and is exactly normalized to unity because of b quantum number conserva- 
tion. It contains the free quark decay as a limiting case with /free(0 = ^{^ ~ f^h/^B)- The 
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distribution function f{^) can be interpreted as the probability of finding a b quark with 
momentum inside the B meson. In addition, we evaluated the mean and variance of the 
distribution function using the techniques of operator product expansion and heavy quark 
effective theory. They specify the primary shape of the distribution function and quantify 
the deviation from the delta function form in the free quark limit. 

The b quark distribution function for the B meson is a key object. Like the well-known 
parton distribution functions for the nucleon in deeply inelastic lepton-nucleon scattering, 
the knowledge of it would help us greatly in understanding the nature of confinement and 
the structure of the B meson. One should try to calculate it using nonperturbative QCD 
methods such as lattice simulation and QCD sum rules. On the other hand, the distribu- 
tion function can be extracted directly from experimental data. The underlying common 
structure of long-distance strong interactions correlates the B Xs'~f and B Xfiu pro- 
cesses. The universality of the distribution function implies that once it is measured from 
one process, it can be used to make predictions in all other processes in a model-independent 
manner. We have discussed the direct extraction of the distribution function from the mea- 
surements of the B — > Xg^ photon spectrum or the B Xfiu spectra dT/d^^. We stress 
that these decay spectra are unique in that they offer the intrinsically most sensitive probe 
of long-distance strong interactions. The experimental extraction of the distribution func- 
tion will lead to a significant improvement of the theoretical description of the B Xsj 
photon energy spectrum, which is very important for seeking new physics. The extracted 
distribution function will also improve the theoretical description of inclusive semileptonic 
B decays, allowing for more precise determinations of \ Vuh\ and \ Vcb\- Moreover, a confronta- 
tion of experimental extraction of the distribution function with QCD predictions will be 
a test of our understanding of the B meson structure and nonperturbative techniques. A 
direct extraction of the distribution function will therefore be an important aspect in future 
experiments in inclusive B meson decays. 

Measurements of 5 — > Xs'~f decays can be used to determine the CKM matrix element 
\Vts\- We have described the theoretically clean methods for the determinations of the 
CKM matrix elements by avoiding the dominant hadronic uncertainties using the known 
normalization of the distribution function or the cancellation of the distribution function 
in the ratio of the differential decay rates. The residual hadronic uncertainty in \Vts\ due 
to higher-twist, power-suppressed corrections is expected to be at the level of one percent. 
With hadronic uncertainties well under control, these methods eventually will yield the 
most accurate value of \Vts\, which is, on the other hand, probably one of the most reliable 
quantities to signal new physics in S — > Xg^. The agreement of \Vts\ extracted from the 
flavour changing neutral current process B — > Xs^ with the value obtained from the direct 
measurements plus unitarity under the assumption that the Standard Model is valid can be 
used to impose constraints on new physics models. 

The universality of the distribution function in radiative and semileptonic inclusive de- 
cays of B mesons is valid only to the leading order in the light-cone expansion. To what 
extent this is a good approximation can be tested experimentally in a variety of ways in the 
B Xg'y, B — > Xuiiy and B XcC-v processes. Higher twist effects may be numerically 
sizable in some region of phase space, especially in the resonance domain. Their quantitative 
consequences deserve a thorough investigation. 

The calculation of long distance effects presented in this paper must be combined with 
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the perturbative QCD and electroweak radiative corrections to give a detailed theoretical 
description of the photon energy spectrum. The theoretical development toward a treatment 
of inclusive radiative decays of B mesons from first principles, in conjunction with precision 
measurements made possible by B factories in the coming few years, will make a decisive 
test of the Standard Model at one-loop level and, more excitingly, might corroborate the 
existence of new physics in inclusive radiative decays of B mesons. 
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